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Abstract 

The J/ip decay angular distributions have been measured in inelastic photopro- 
duction in ep collisions with the ZEUS detector at HERA, using an integrated 
luminosity of 468 pb~^. The range in photon-proton centre-of-mass energy, W , 
was 50 < IV < 180 GeV. The J /if) mesons were identified through their de- 
cay into muon pairs. The polar and azimuthal angles of the ix^ were measured 
in the J/ip rest frame and compared to theoretical predictions at leading and 
next-to- leading order in QCD. 
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1 Introduction 



In the HERA photoproduction regime, where the virtuahty of the exchanged photon is 
small, the production of inelastic J/il) mesons is dominated by boson-gluon fusion: a 
photon emitted from the incoming lepton interacts with a gluon coming from the pro- 
ton to produce a cc pair which subsequently forms & J/ip meson. Production of J/ip 
through boson-gluon fusion can be calculated using perturbative Quantum Chromody- 
namics (pQCD) in the colour-singlet (CS), in the non-relativistic QCD (NRQCD) [1,2] 
and in the /cT-factorisation frameworks [3,4]. 

In the CS approach, only the colourless cc pair produced in the hard subprocess can lead 
to a physical J/il) state. In the NRQCD approach, a cc pair emerging from the hard 
process in a colour-octet (CO) state can also evolve into a, J/ip state with a probability 
proportional to universal long-distance matrix elements (LDME) that are obtained from 
experiment. In the /cT-factorisation approach, the effects of non-zero incoming parton 
transverse momentum are taken into account. Cross sections are then calculated in the 
CS approach as a convolution of unintegrated (transverse-momentum dependent) parton 
densities and leading-order (LO) off-shell matrix elements. 

At HERA, measurements of inelastic J /ip differential cross sections [5,6] are reproduced by 
a next-to-leading-order (NLO) QCD calculation [7,8] performed in the CS framework. The 
measurements are also reasonably well described by LO CS plus CO calculations [9-11], 
with LDME as determined in a LO analysis of hadroproduction and 5-decay data [12-15]. 

The polar and azimuthal distributions of the J/il) decay leptons in the J/il) rest frame 
may be used to distinguish between CS and CO models. These helicity distributions 
are expected to be different as a function of the J/ip transverse momentum, px-, and 
inelasticity, the fraction of the incident photon energy carried by the J/ if) in the proton 
rest frame, z [16]. 

Helicity-distribution measurements have already been performed by the ZEUS [5] and 
HI [6] collaborations. In the final study presented here, J /ip mesons were identified using 
the decay mode J /ip f^~^f^~~ and were measured in the range 50 < < 180 GeV, where 
W is the photon-proton centre-of-mass energy. The data sample under study includes the 
data used in the previously published ZEUS analysis [5] and corresponds to an increase 
in statistics of a factor of 12. 
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2 Experimental set-up 



The analysis presented here is based on data collected by the ZEUS detector at HERA 
in the period 1996-2007. In 1998-2007 (1996-1997), HERA provided electroif] beams of 
energy = 27.5 GeV and proton beams of energy Ep = 920 (820) GeV, resulting in a 
centre-of-mass energy of y/s = 318 (300) GeV, corresponding to an integrated luminosity 
of 430 ± 11 (38 ± 0.6) pb-i. 

A detailed description of the ZEUS detector can be found elsewhere [17,18]. A brief 
outline of the components that are most relevant for this analysis is given below. Charged 
particles were tracked in the central tracking detector (CTD) [19], which operated in 
a magnetic field of 1.43 T provided by a thin superconducting coil. Before the 2003- 
2007 running period, the ZEUS tracking system was upgraded with a silicon microvertex 
detector (MVD) [20]. The high-resolution uranium-scintillator calorimeter (GAL) [21] 
consisted of three parts: the forward (FGAL), the barrel (BGAL) and the rear (RGAL) 
calorimeter^. 

Muons were identified by tracks in the barrel and rear muon chambers (BMUON and 
RMUON) [22]. The muon chambers were placed inside and outside the magnetised iron 
yoke surrounding the GAL. The barrel and rear inner muon chambers (BMUI and RMUI) 
covered the polar-angle regions 34° < 6 < 135° and 135° < 9 < 171°, respectively. The 
luminosity was measured using the Bethe-Heitler reaction ep e'jp with the luminos- 
ity detector which consisted of a lead-scintillator calorimeter [23] and, after 2002, an 
additional magnetic spectrometer [24] system. 



3 Event selection 

Inelastic events are often selected using the inelasticity, z. In this analysis, however, the 
events were selected using the transverse momentum, pt, of the J/ip and additional ac- 
tivity in the detector. This kind of selection permits direct comparisons with the different 
theoretical predictions [16]. 

The online and offline selections as well as the reconstruction of the kinematic variables 
closely follow a previous analysis [5]. 

^ Here and in the following, the term "electron" denotes generically both the electron (e^) and the 
positron (e"*"). 

^ The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" . and the X axis pointing towards the 
centre of HERA. The coordinate origin is at the nominal interaction point. The polar angle, 9, is 
measured with respect to the proton beam direction. The pscudorapidity is defined as r;=-ln(tan |). 
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Online, the BMUI and RMUI chambers were used to tag muons by matching segments 
in the muon chambers with tracks in the CTD/MVD, as well as with energy deposits in 
the GAL consistent with the passage of a minimum- ionising particle (m.i.p.). 

Offline, an event was accepted if it had two tracks forming a. J /ip candidate. One track had 
to be identified in the inner muon chambers and matched to a m.i.p. cluster in the GAL. 
It was required to have a momentum greater than 1.8 GeV if it was in the rear region, 
or a transverse momentum greater than 1.4 GeV if in the barrel region. The other track 
had to be matched to a m.i.p. cluster in the GAL and was required to have a transverse 
momentum greater than 0.9 GeV. Both tracks were restricted to the pseudorapidity region 
|?7| < 1.75. To reject cosmic rays, events in which the angle between the two muon tracks 
was larger than 174° were removed. 

The pt of the J / ip candidate was required to be larger than 1 GeV. In addition, events were 
required to have an energy deposit larger than 1 GeV in a cone of 35° around the forward 
direction (excluding possible calorimeter deposits due to the decay muons). According 
to Monte Garlo (MG) simulations, these requirements completely reject exclusively pro- 
duced J/ip mesons {ep — > epJ/ip) as well as proton-diffractive events {ep — > eYJ/%lj) in 
which the mass of the proton dissociative state. My, is below 4.4 GeV. To further reduce 
diffractive background, events were also required to have at least one additional track 
with a transverse momentum larger than 0.125 GeV and pseudorapidity |?7| < 1.75. 



4 Kinematic variables and signal extraction 

The photon-proton centre-of-mass energy, W , is: 

W^={P + q)\ (1) 

where P and q are the four-momenta of the incoming proton and exchanged photon, 
respectively. It was calculated using: 

W^ = 2E,{E-pz) (2) 

where {E — pz) is summed over all final-state energy-flow objects [25] (EFOs) which 
combine the information from calorimetry and tracking. 

The inelasticity z = was determined as: 

_ {E-pz)j/^ 

where pj/^ is the four-momentum of the J/ip and {E — pz)j/^ was calculated using the 
tracks forming the J /if). 
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The kinematic region considered was 50 < 14^ < 180 GeV where the acceptance was 
always above 10%. A requirement of E — pz < 20 GeV restricted the virtuahty of the 
exchanged photon = —q^ < 1 GeV^, with a median of ~ 10^"^ GeV^. The ehmination 
of deep inelastic scattering events was independently confirmed by searching for scattered 
electrons in the GAL [26]; none was found. 

The invariant-mass spectrum of the J/ip candidates with > 1 GeV and 2 > 0.1 is 
shown in Fig.[Tl Both the J /if) and if)' peaks are visible. The background was estimated by 
fitting the product of a second-order polynomial and an exponential to the region outside 
the invariant-mass window, 2.85-3.3 GeV. The number of J/ip events was obtained by 
subtracting the number of background events estimated from the fit procedure from the 
total number of events inside the invariant-mass window; 12310 ± 140 J /if) events were 
found. As the signal to background ratio is large, the extracted number of J/ip events 
has little sensitivity to the analytical form of the function used for the background fit. 

5 Monte Carlo and background evaluation 

The inelastic production of J /if) mesons was simulated using the Herwig 5.8 [27] MG 
generator, which generates events according to the LO diagrams of the boson-gluon- 
fusion process, 7(7 — J/ipg, as calculated in the framework of the GS model. This process 
is called a direct photon process, because the incoming photon couples to the c quark 
directly. The Herwig MG sample was reweighted in pt, z and W in order to give the 
best description of the data. 

There are other sources of J/ "dp mesons which were classified as background in the present 
analysis and were estimated either from MG models or previous measurements. Although 
the relative rate of each process is given below, the helicity distributions of these Jji\} 
sources are poorly known, so the contributions were not subtracted. 

Diffractive production of Jj^l) mesons with proton dissociation was simulated with the 
Epsoft [28] MG generator, which was tuned to describe such processes at HERA [29]. 
This background is suppressed by the requirement on the tracks and by the cut on the 
minimum pt of the J/^/'. The overall contribution of this background is 6%; it is largest 
in the lowest pt bin (1 < Pt < 1-4 GeV), where it is 7.5%, and in the highest z bin 
(0.9 < z < 1), where it is 66%. 

The production of J/ ^/^ mesons originating from B-meson decays was simulated using the 
Pythia 6.2 MG generator [30]. The beauty-quark mass was set to 4.75 GeV and the B 
to J l^l) branching ratio was set to the PDG value [31]. According to the MG, 1.6% of the 
observed J/?/' events were from 5-meson decays; the fraction is largest in the highest pj^ 
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bin (4.2 < < 10 GeV), where it is equal to 6.3%, and in the lowest z bin (0.1 < z < 0.4), 
where it is 8.4%. 

The background from ^Z^' to J j'^l) decays is expected to be around 15%, as obtained using 
the direct measurement of the V^' to J/ 1^) cross section ratio [5] and the branching ratio of 
the to J jil). 

All generated events were passed through a full simulation of the ZEUS detector based on 
Geant 3 [32]. They were then subjected to the same trigger requirements and processed 
by the same reconstruction program as the data. 



6 Reconstruction of the helicity parameters 

The helicity analysis was performed in the so-called "target frame" [16], i.e. the J/V' rest 
frame with the axes Z' = —Z and Y' along the vector q' x [—P'), where q' and P' are 
the three-vectors associated with the exchanged photon and incoming proton. The polar 
and azimuthal angles of the /i+ in this frame are denoted 6* and (p*. 

The differential cross sections in 6* and can be parametrised as [16]: 

da 
d cos 9* 

and 



ocl + Acos'r, (4) 



da X u 

_«l + _ + _„,2r, (5) 

where A and u, the polar and azimuthal angular parameters, are functions of px and z. 
The predictions for A and u depend on the production mechanism. The value A = +1 
corresponds to J/tp mesons fully transversally polarised, while A = — 1 corresponds to 
J/ip mesons fully longitudinally polarised. 

The A and u parameters were determined in bins of z and pt, each time integrating over 
the other variable. As a function of pt, the integration range for z was set to 0.4 < z < 1, 
thereby avoiding the region 0.1 < z < 0.4 where the ratio of signal to combinatorial 
background is rather poor (0.52). The integration range in px started at = 1 GeV. 

In the estimation of the parameters A and u, the helicity distributions of the background 
events present under the J/ip peak were added to the MC distributions. The shape of 
the background helicity distributions was taken from the side bands, while the number of 
background events was taken from the fits described in Section HI 

The Herwig MG generator-level distributions dN/d cos 6* {dN/d(j)*) were re-weighted 
according to Eq. |4] ([5]) within a search grid of A (z/) values. For each re-weighted distri- 
bution, the value of '^as calculated from a comparison to the data. The A (z/) value 
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providing the minimum x^i Xmin' taken as the central value. The parameter values 
with = Xmin + 1 were used to calculate the statistical uncertainties. The Xmin P^r 
degree of freedom were typically around one. Equation H] was first used to extract A, and 
then A was inserted into Eq. [5] to extract u, see Tables [H and O 

7 Systematic uncertainties 

The following sources of systematic uncertainties were investigated (their effects are given 
in parentheses): 

• muon chamber efficiencies: the BMUI and RMUI muon chamber efficiencies were 
extracted from the data using muon pairs coming from elastic J/ip events and from 
the process 77 fi'^fi". These efficiencies are known up to an uncertainty of about 
±5% (< 5% of the statistical error); 

• analysis cuts: this class comprises the systematic uncertainties due to the uncertainties 
in the measurement of momentum, transverse momentum and pseudorapidity of the 
muon decay tracks. Each cut was varied within a range determined by the resolution 
in the appropriate variable (< 5% of the statistical error); 

• CAL energy scale: the CAL energies were varied by ±5% in the simulation, in accor- 
dance with the uncertainty in the CAL energy scale (on average 10% of the statistical 
error) ; 

• hadronic energy resolution: the W and z resolutions are dominated by the hadronic 
energy resolution affecting the quantity {E — pz)- The [E —pz) resolution in the MC 
was smeared event by event by ±20% (on average 10% of the statistical error); 

• pt, W and z spectra: the pt, W and z spectra of the J/z/^ mesons in the Herwig 
MC simulation were varied within ranges allowed by the comparison between data and 
simulation (on average 15% of the statistical error); 

• additional track requirement: the kinematic cuts for the additional track requirement 
were tightened and loosened in both data and MC (on average 15% of the statistical 
error) ; 

• influence of diffractive contamination at high z on the A and i> extractions as a function 
oi Pt'- a and u were extracted changing the z integration range from 0.4 < 2; < 1 to 
0.4 < z < 0.9 (on average 25% of the statistical error); 

• angular coverage in 6*: cos 6** was restricted to the range —0.8 < cos^* < 0.8 in order 
to avoid low-acceptance regions (on average 30% of the statistical error); 
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• invariant-mass window: the J / ip invariant mass window was enlarged by 50 MeV and 
tightened by 100 MeV (on average 30% of the statistical error). 

All of the above individual sources of systematic uncertainty were added in quadrature. 
No systematic uncertainties are quoted for J/ip coming from 5-meson decays and J/ if) 
coming from if)' decays. The uncertainties on the integrated-luminosity determination 
and on the J/if) ^ f^^f^~ branching ratio, which would result in an overall shift of a cross 
section measurement, do not contribute to the measurement of the helicity parameters. 



8 Results 

The values of the parameter A are shown as a function of pt and z in Fig. [2ti) and b), 
respectively; the values of the parameter z/ are displayed in Fig. [3^) and b). All the values 
are also listed in Tables [1] and [2l The data indicate that the parameter A depends, if at 
all, only weakly on px and rises slowly with z. The parameter u does not seem to depend 
on Pt, while it seems to increase at low and high z. 

The data are compared to various theoretical predictions for photoproduction at Q'^ = 0. 
These predictions do not consider the polarisation due to J/ 4' coming from ijj' decays, 
5-meson decays and from diffractive processes. The curves identified by the label LO CS 
show the LO prediction in the CS framework including both direct and resolvecjfl processes. 
The two lines identified by the label LO + /ct (JB) and LO + /ct (dGRV), represent 
the predictions of a /cj^-factorisation model [33] using two different unintegrated gluon 
distributions and including only direct processes. The band identified by the label NLO 
CS represents the predictions of a NLO calculation [34,35] including only direct processes. 
The width of the band gives the uncertainties of the calculation due to variations of the 
renormalisation and factorisation scales. It stops aX z = 0.9 because no reliable predictions 
can be obtained near 2; = 1 for this fixed-order calculation. The band, identified by 
the label LO CS+CO, shows the LO prediction [16] including both CS and CO terms, 
including both direct and resolved processes. The width of the band results from the 
uncertainties in the values of the long-distance matrix elements. In Fig. [2^) and [3^), with 
z integrated up to 2; = 1, the LO CS+CO cross section is CO dominated. 

None of the models provides predictions for both A and v that agree well with the data 
everywhere in z and pt- For A as a function of pt-, all the models roughly describe the 
data, with the NLO CS prediction providing the poorest description. For high values 
of Pt, the polarisation in the data remains small, while LO CS predicts a progressive 

In these resolved processes, the incoming photon does not couple to the c quark directly, but via its 
hadronic component. They are expected to contribute mainly to the region of z < 0.4. 
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increase and NLO CS and LO a progressive decrease. The LO CS + CO prediction 
remains flat, with a small and positive value of A. For A as a function of z, all theoretical 
predictions are in rough agreement with the data. The pt and z dependencies of u are not 
described by the LO CS predictions, while the other models provide better descriptions 
of the data. 

The NLO CS calculation for pt > I GeV suffers from large scale uncertainties connected 
to the presence of negative values of the diagonal components of the spin density matrix 
Pt ^ ^ GeV [34,35]. In order to avoid this problem, measurements and calculations 
were repeated increasing the pt cut flrst to 2 GeV and then to 3 GeV. In Fig. H^) and 
b) the A and v parameters, respectively, are shown as a function of z for px > 2 GeV, 
while in Fig. Ht) and d) the same parameters are displayed for p^- > 3 GeV. All the 
values are listed in Table [3l The NLO CS calculation [35], also shown in these figures, has 
now smaller uncertainties, but the agreement with the data is only satisfactory for the i' 
parameter. Sizeable discrepancies remain for the A parameter both for pt > 2 GeV and 
pr > 3 GeV. 

9 Conclusions 

The J/ ip helicity distributions in the inelastic photoproduction regime have been measured 
using a luminosity of 468 pb~^. The J/ip helicity parameters A and u were extracted in 
the target frame as a function of the transverse momentum and of the inelasticity of the 
J/ip. The results were compared to LO QCD predictions in the colour-singlet, colour- 
singlet plus colour-octet and factorisation frameworks. A recent NLO QCD prediction 
in the colour- singlet framework was also considered. Even though the experimental and 
theoretical uncertainties are large, none of the predictions can describe all aspects of the 
data. 
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Pt range (GeV) 



1.0-1.4 
1.4-1.9 
1.9-2.4 
2.4-3.4 
3.4-4.2 
4.2-10. 



(pt) (GeV) 



1.2 
1.6 
2.1 
2.8 
3.7 
5.2 



A 



0.161°: 

0.38l[J: 

-0.151°: 
0.211°: 

0.341° 



-0.10 
-0.04 



0.31 



17 +0.10 
17 -0.10 
0.20 
0.08 
0.26 
0.05 
0.16 
0.17 
0.21 
0.19 



o.26i[!: 
-0.121°: 
-o.33lg: 
-0.09115: 



0.21: 

-0.50 



10 +0.05 
10 -0.08 
.12 +0.07 
.14 -0.01 
.21 +0.12 
.24 -0.16 
.19 +0.13 
.20 -0.13 
26 +0.10 
28 -0.03 
,27 +0.11 
.28 -0.07 



Table 1: J/ip helicity parameters A and u as a function of pt measured in the 
target frame for 50 < W < 180 GeV, O.4 < z < 1 and px > I GeV. The first 
uncertainty is statistical and the second is systematic. 



z range 


{^) 


A 


V 


0.10-0.40 
0.40-0.55 
0.55-0.70 
0.70-0.80 
0.80-0.90 
0.90-1.0 


0.27 
0.48 
0.61 
0.75 
0.85 
0.95 


n 97+0.43 +0.13 
^■^ ' -0.40 -0.14 
n 99+0.19 +0.18 
'-'•'^^-0.18 -0.06 
n -[C+0.12 +0.10 
"-'•-'-"-0.13 -0.04 
n qQ+0.15 +0.04 
"-'•'-'^-0.15 -0.13 

n Qn+0-17 +0.15 

'-'■'-"-'-0.17 -0.05 
n /IQ+0.39 +0.15 
'-'•^'^-0.34 -0.08 


n QQ+0-24 +0.14 
U.yj_Q 27 _o,40 

n. m +0.17 +0.08 
'-'•'-'-'--0.18 -0.03 

n nc;+o-io +0.05 
'-'•'-'^-0.11 -0.02 

n. 07+0.11 +0.04 
'-'•'-"-0.13 -0.07 

n n4+o-i4 +0.08 

'-'•'-'^-0.15 -0.08 
n cr/i+0.26 +0.16 
'-'•^^-0.28 -0.08 



Table 2: J helicity parameters A and v as a function of z measured in the 
target frame for 50 < W < 180 GeV, 0.1 < z < 1 and pr > 1 GeV. The first 
uncertainty is statistical and the second is systematic. 



Pt range (GeV) 


z range 


(z) 


A 


U 


Pt>2 


0.10-0.55 
0.55-0.70 
0.70-0.80 
0.80-0.90 


0.37 
0.61 
0.75 
0.85 


n qp;+0.34 +0.04 
'-'•'-'^-0.31 -0.41 

n nc;+o-i8 +0.20 

'-'•'-'^-0.17 -0.06 
n q/i+0.22 +0.04 
'-'•'-'^-0.22 -0.17 
n 1 9+0.24 +0.29 
'-'•1-^-0.23 -0.08 


n 07+0.35 +0.08 
'-'•'-" -0.39 -0.09 
n 1 r,+0.19 +0.08 
'-'•-'-'-'-0.20 -0.05 

n nc;+o-2i +0.06 

U.UO_Q23 _o.25 

n qQ+0.26 +0.22 
u.oy„Q 27 „o.01 


Pt > S 


0.10-0.55 
0.55-0.70 
0.70-0.90 


0.38 
0.62 
0.79 


n 8n+0-53 +0.11 
U.OU_o.45 _o.43 

n 9(3+0.31 +0.23 
'-'•^'-'-0.28 -0.06 

n nQ+0-25 +0.27 
U.UJ_g 23 _o.l2 


n (-,7+0-40 +0.07 
^■'J 1 -0.44 -0.04 

n 9(3+0.27 +0.20 
'-'•^'-'-0.28 -0.06 

n qp:+0.27 +0.09 
U.OJ„Q 28 „o.07 



Table 3: J/tp helicity parameters A and u as a function of z measured in the 
target frame for 50 < W < 180 GeV, 0.1 < z < 1 and for pt > 2 and pt > 3 GeV. 
The first uncertainty is statistical and the second is systematic. 
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Figure 1: Dimuon invariant mass, m^^, spectrum in the phase-space region 
50 < < 180 GeV, z > 0.1 and pt > 1 GeV. The continuous line represents the 
fitted background. 
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Figure 2: The helicity parameter X, measured in the target frame, as a function 
of (a) Pt, and (b) z. The measurement is performed in the kinematic range 50 < 
W < 180 GeV, 0.1 < z < 1 and px > 1 GeV. The measurement as a function 
of Pt is restricted to the kinematic range 0.4 < z < 1. The inner (outer) error 
bars correspond to the statistical (total) uncertainty. The theoretical curves are 
described in the text. 
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Figure 3: The helicity parameter v, measured in the target frame, as a function 
of (a) Pt, o,nd (b) z. The measurement is performed in the kinematic range 50 < 
W < 180 GeV, 0.1 < z < 1 and px > 1 GeV. The measurement as a function 
of Pt is restricted to the kinematic range 0.4 < z < 1. The inner (outer) error 
bars correspond to the statistical (total) uncertainty. The theoretical curves are 
described in the text. 
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Figure 4: Distributions of the helicity parameters (a), (c) X and (b), (d) v as a 
function of z, measured in the target frame, for 50 < TV < 180 GeV, 0.1 < 2; < 0.9 
and (a), (b) pt > 2 GeV and (c), (d) px > 3 GeV. The inner (outer) error bars 
correspond to the statistical (total) uncertainty. The theoretical bands are described 
in the text. 
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